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ABSTRACT 

This paper concerns the K band Hubble diagram for the brightest cluster galaxies 
(BCGs) in a sample of X-ray clusters covering the redshift range 0.05 < z < 0.8. We 
show that BCGs in clusters of high X-ray luminosity are excellent standard candles: 
the intrinsic dispersion in the raw K band absolute magnitudes of BCGs in clusters 
with Lx > 2.3 x 10 44 erg s _1 (in the 0.3 - 3.5 keV band) is no more than 0.22 mag, 
and is not significantly reduced by correcting for the BCG structure parameter, a, 
or for X-ray luminosity. This is the smallest scatter in the absolute magnitudes of 
any single class of galaxy and demonstrates the homogeneity of BCGs in high-Lx 
clusters. By contrast, we find that the brightest members of low-Lx systems display a 
wider dispersion (~ 0.5 mag) in absolute magnitude than commonly seen in previous 
studies, which arises from the inclusion, in X-ray flux-limited samples, of poor clusters 
and groups which are usually omitted from low redshift studies of BCGs in optically 
rich clusters. 

Spectral synthesis models reveal the insensitivity of K band light to galaxy evolu- 
tion, and this insensitivity, coupled with the tightness of its Hubble relation, and the 
lack of evidence of significant growth by merging (shown by the absence of a correla- 
tion between BCG structure parameter, a, and redshift), makes our sample of BCGs 
in high-Lx clusters ideal for estimating the cosmological parameters f^M and f^A, free 
from many of the problems that have bedevilled previous attempts using BCGs. The 
BCGs in our high-Lx clusters yield a value of 57m = 0.28 ± 0.24 if the cosmological 
constant A = 0. For a flat Universe we find f^M = 0-55io'i5 with a 95 per cent con- 
fidence upper limit to the cosmological constant corresponding to f^A < 0.73. These 
results are discussed in the context of other methods used to constrain the density of 
the Universe, such as Type la supernovae. 
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1 INTRODUCTION 

The Hubble (redsfiift-magnitude) diagram for brightest 
cluster galaxies (BCGs) is a classic cosmological tool (re- 
cently reviewed by Sandage 1995). Sandage and collabora- 
tors (e.g. Sandage 1972a,b; Sandage & Hardy 1973; Sandage, 
Kristian & Westphal 1976) used the Hubble diagram for 
BCGs to verify Hubble's linear redshift-distance law out 
to a distance one hundred times greater than that probed 
by the initial bright galaxy sample of Hubble (1929), and 
attempted to measure the deceleration parameter from its 
deviation from the linear Hubble law. More recently, the 
BCG redshift-magnitude relation has appeared in studies of 
the formation and evolution of elliptical galaxies (Aragon- 



Salamanca et al. 1993) and of large-scale bulk flows in the 
Universe (Lauer & Postman 1994; Postman & Lauer 1995; 
Hudson & Ebeling 1997). 

Underlying all this work is the fact that BCGs appear 
to be good standard candles after correction for systematic 
effects and the dependence of absolute magnitude on galaxy 
structure and environment, with their aperture luminosi- 
ties at low redshift showing a dispersion of only ~ 0.3 mag 
(Sandage 1988). This small dispersion is all the more re- 
markable since BCGs do not appear to be a particularly 
homogeneous class of objects at first sight: some exhibit the 
large, extended envelopes that make them cD galaxies; many 
(up to half according to Hoessel & Schneider 1985) have mul- 
tiple nuclei; and their immediate local environments vary 
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from relative isolation in the cluster core to interaction and 
merging, in the case of dumbbell galaxies. BCGs tend to 
dominate their host clusters visually, but little is known of 
their place in the formation and evolution of the clusters, or 
of their own origins. 

Answers to some of these questions may result from ex- 
tending the Hubble diagram for BCGs to higher redshifts, 
but there are several problems with that. Firstly, as one 
moves to higher redshifts, optical wavebands become in- 
creasingly sensitive to the effects of star formation, and K- 
corrections become large and uncertain. This problem may 
be circumvented, as discussed in detail below, by working 
in the K band, which continues principally to sample the 
mature stellar population out to high redshifts: the K band 
Hubble diagram for a small sample of BCGs was presented 
by Aragon-Salamanca et al. (1993). 

A more serious problem is that of cluster selection. In 
the past, cluster samples have been selected from optical 
surveys, on the basis of identifying enhancements in the 
surface density of galaxies above a fluctuating background. 
Such a procedure suffers from the well-known projection 
effects caused by chance alignments of galaxies along the 
line-of-sight (van Haarlem, Frenk & White 1997, and refer- 
ences therein). This contamination will be more severe at 
fainter magnitudes, thus strongly discouraging the use of 
this method for detecting high-redshift clusters. A solution 
to this problem is provided by X-ray cluster selection. Most 
bona fide clusters contain hot gas in their cores, emitting 
a copious flux of X-rays and making them visible to large 
distances. The intensity of the thermal bremsstrahlung X- 
ray emission from a cluster is directly related to the depth 
of its gravitational potential well, unlike its optical galaxy 
richness, while the compactness of the X-ray-emitting re- 
gion, compared to the extent of the galaxy distribution, also 
means that projection effects are minimal in comparison to 
those arising in optical cluster selection (Romer et al. 1994; 
van Haarlem et al. 1997). 

These concerns motivated the current study, the long- 
term goal of which is to investigate the formation and evolu- 
tion of BCGs in host clusters with known X-ray luminosities, 
by studying their physical properties over a wide span of cos- 
mological time. In this paper we concentrate on extending 
the infrared Hubble diagram for BCGs to higher redshifts 
than studied before using K band images of a large sample of 
BCGs in X-ray clusters with z ^ 0.8. The plan of the remain- 
der of this paper is as follows. In Section 2 we briefly describe 
our cluster sample, the observations made and their reduc- 
tion, leading to the presentation of the uncorrected K — z 
diagram for our sample in Section 3, followed by a discus- 
sion of K-correction models for our BCGs in Section 4. We 
consider the physical properties of the BCGs and their host 
clusters in Section 5 and look for correlations between them 
in Section 6. These correlations are used to produce a cor- 
rected Hubble diagram in Section 7, which is used to derive 
constraints on the cosmological parameters Qm and J1a in 
Section 8. In Section 9 we discuss the results of this paper 
and present the conclusions we draw from them. An Ap- 
pendix tabulates some statistical results omitted from the 
main body of the text. 
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Figure 1. The luminosities of the EMSS clusters (in the 0.3-3.5 
keV passband) against redshift. The filled circles are the clusters 
imaged at K and discussed in this paper, while the empty cir- 
cles are the rest of the EMSS cluster sample of Gioia & Luppino 
(1994). The two triangles are the ROSAT clusters. Note that 
there is no single, sharp flux limit to the EMSS survey — clus- 
ter detections were based on pointed observations whose limiting 
sensitivity ranges from 5 X 10~ 14 erg cm -2 s" 1 to 3 X 10 — 12 erg 



2 THE DATA 

2.1 The X-ray Cluster Sample 

The X-ray clusters for this work were chosen principally from 
the Einstein Medium Sensitivity Survey (Gioia et al. 1990, 
Stocke et al. 1991, hereafter EMSS) catalogue of 104 X-ray 
selected clusters (Gioia & Luppino 1994), but supplemented 
by two clusters from the ROSAT All-Sky Survey (Voges 
1992, Triimper 1993). Fig. Ill shows the X-ray luminosities 
listed in Gioia & Luppino (1994) for all the EMSS clusters, 
as a function of redshift, compared with the sample studied 
here: the X-ray luminosities were calculated assuming that 
Qm = 1, Qa — and the Hubble constant is Ho = 50 km 
s _1 Mpc -1 . No attempt was made to observe a statistically 
complete subsample of EMSS clusters, although it is clear 
from Fig. hi that those selected are broadly representative of 
the range of X-ray luminosities for clusters within the parent 
survey: a two-sided Kolmogorov-Smirnov (K-S) test yields 
a probability of 0.26 that the cumulative X-ray luminosity 
distributions of the EMSS clusters imaged and not imaged 
would differ by more than observed if both subsamples were 
drawn from the same parent distribution. The corresponding 
K-S test on the redshift distributions of EMSS cluster sub- 
samples imaged and not imaged yields a probability of only 
0.04, reflecting the fact that we preferentially selected EMSS 
clusters at higher redshifts, where the luminosity distance 
varies appreciably with cosmology and, thus, constraints on 
Qm and Ha may be derived from the BCG Hubble diagram. 



2.2 Observations 

Observations were made on 8-10 November 1994 and 20- 
22 April 1995 using the IRCAM3 infrared camera on the 
3.8m United Kingdom Infrared Telescope (UKIRT): IR- 
CAM3 contains a 256 x 256 InSb array, with a pixel size 
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of 0.286 arcsec. Two observing procedures were followed. 
For distant clusters, a series of 100s (10 x 10s co-added) 
exposures were taken, with the galaxy offset in a five-point 
jitter pattern within the initial field of view of the chip. For 
nearby clusters, where the BCG covered too large an area 
of the chip to allow accurate sky level determination by this 
method, sky exposures were obtained after each galaxy ex- 
posure, allowing the jitter pattern for the sky frames (four-, 
or eight-point) to be built-up simultaneously with that of 
the jitter pattern for the galaxy frames (five-, or nine-point). 
Typical total on-source integration times used were ~ 1000s 
for z <, 0.2, 1500-3000s for 0.2 <, z <> 0.4 and 4500s above 
z — 0.5. Dark frames were taken regularly throughout the 
nights, as were standard star observations, using sources 
from the UKIRT Faint Standards list (Casali & Hawarden 
1992). 
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2.3 Data Reduction 

The data were reduced using standard iraft] routines and 
the following procedure. Appropriate dark frames were sub- 
tracted from each galaxy or sky image, and the resulting 
frames were multiplied by a mask, flagging known bad pixels 
in the array to be ignored subsequently, before normalising 
the frames to unit median: we found the median to be a 
more stable measure of the sky level than the mode. Flat- 
field frames were constructed, by median filtering the galaxy 
frames (for distant clusters) or sky frames (for nearby clus- 
ters), and these flat frames were normalised to unit median. 
The galaxy frames were then divided by the appropriate 
flat field frame and final mosaic images produced by cross- 
correlating the individual flat-fielded galaxy frames to de- 
termine their relative offsets. A similar procedure was used 
to reduce the standard star frames. 

A total of 52 BCGs were observed. Two of these, (those 
in MS1401.9+0437 and MS1520.1+3002), did not yield pho- 
tometric quality data and were dropped from our sample, 
as were two others (MS1209.0+3917 and MS1333.3+1725) 
which have dubious redshift identifications in Gioia & Lup- 
pino (1994). A fifth BCG, in MS0440. 5+0204, was sub- 
sequently omitted, because it is close to a bright star, 
which falls within our chosen photometric aperture. This 
left us with a sample of 47 BCGs: note that Stocke et al. 
(1991) and Gioia & Luppino (1994) suggest that the two 
EMSS sources MS2215. 7-0404 and MS2216.0-0401 may be 
part of the extended X-ray emission from a single cluster, 
but we have treated them as two distinct objects. We as- 
sume cosmological matter and vacuum energy densities of 
$!m = 1 and £7a = 0, respectively, and h = 0.5 (where 
Hq = 100ft kms~ 1 Mpc~ 1 is the Hubble constant) through- 
out this paper, unless stated otherwise. 



3 THE RAW K-z DIAGRAM 

We obtained K band aperture magnitudes for our BCGs 
using the iraf package APPHOT, with a fixed metric aperture 



IRAF is distributed by National Optical Astronomy Observa- 
tories, which is operated by the Association of Universities for 
Research in Astronomy, Inc., under contract to the National Sci- 
ence Foundation. 



Figure 2. The raw K — z diagram. The filled circles show the 
apparent K band magnitudes of the 47 BCGs in our X-ray sample, 
while the empty circles are the BCGs from the optical sample 
of Aragon-Salamanca et al. (1993), excluding the one BCG in 
common between the two samples. The magnitudes are corrected 
for Galactic extinction, but have not been K-corrected. 



of diameter 50 kpc: in the few cases, where there were two 
similarly bright galaxies in the cluster, the BCG was taken 
to be the brighter of the two as judged by this magnitude. 
Airmass corrections were made using the standard median 
Mauna Kea K band extinction value of 0.088 mag/airmass. 

The choice of aperture was made on the basis of two 
considerations. Firstly, it allows direct comparison with the 
results of Aragon-Salamanca et al. (1993) who studied BCGs 
in nineteen optically-selected clusters. Secondly, as discussed 
in Section 5.1, we require a sufficiently large metric aperture 
to allow accurate estimation of the BCG structure parame- 
ter, a, for our more distant clusters, given the seeing condi- 
tions obtaining during our runs. 

From comparing frames of the same object taken at 
different times on the same night, on different nights of the 
same run and between the two runs, we obtain formal sta- 
tistical errors of 0.04 to 0.07 mags on our 50 kpc magni- 
tudes, except for 11 BCGs observed on one particular night 
for which greater uncertainty in our photometric calibra- 
tion produced uncertainties of 0.1 mags. For the one object 
(MS0015. 9+1609 = C10016+16) we have in common with 
Aragon-Salamanca et al. (1993) they measure an aperture 
magnitude that agrees with ours to well within these es- 
timated errors (15.56 mag compared to our 15.58 mag). In 
addition to these statistical errors, there are likely to be sys- 
tematic errors resulting, for example, from contamination 
from other cluster galaxies falling within the photometric 
aperture: we make no attempt to estimate these (necessarily 
very uncertain) corrections. The aperture magnitude data 
set is given in Table 2 and the resultant raw K — z diagram 
is shown in Fig. H It is clear from Fig. that there is a large 
scatter in the raw K — z diagram for the BCGs in our X-ray 
cluster sample, and that most of the scatter comes from the 
lower- z half of the sample: before expressing that scatter in 
terms of the rms dispersion in the absolute magnitudes of 
the galaxies, and investigating its origin, we must discuss 
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how to K-correct the galaxy magnitudes, which we do in 
the next Section. 



4 INTERPRETING THE BCG K BAND 
MAGNITUDES 

There is a consensus developing that cluster ellipticals are 
old, for example: Chariot & Silk (1994) model the evolu- 
tion of spectral indicators of E/SO galaxies in low-redshift 
(z < 0.4) clusters, and suggest that a couple of per cent at 
most of their stars have been made in the past 2.5 Gyr; 
Bender, Ziegler & Bruzual (1996) use Mgt — a data for 
16 cluster ellipticals in three clusters at z ~ 0.37 to show 
that the majority of their stars must have formed at z > 2, 
and that the most luminous galaxies may have formed at 
z > 4; and Ellis et al. (1997) use rest-frame UV-optical 
photometry of ellipticals in z ~ 0.5 clusters obtained with 
the Hubble Space Telescope to show that the bulk of the 
star formation in the ellipticals in dense clusters was com- 
pleted before z ~ 3 in conventional cosmologies. If we as- 
sume that the same goes for our BCGs as for cluster el- 
lipticals in general, then our photometric data are easy to 
interpret, since, by choosing the K band, we are primarily 
sampling the passively evolving mature stellar population 
over our full redshift range, and the resultant K-correction 
will be insensitive to the exact age of the galaxy: we take the 
K-correction, K(z), as encompassing the bandwidth, band- 
shifting and evolution terms (e.g. Sandage 1995), defining it 
so that the absolute magnitude, M, of a galaxy observed to 
have an apparent magnitude m at redshift z is given by 
M = m + 5 — 51og 10 (dL) — K (z), where dt,(z,riM,£lA, Ho) 
is the luminosity distance (in pc), and where we have ne- 
glected the effect of Galactic extinction. 

In the light of these results, we take as our default K- 
correction model that resulting from a Bruzual & Chariot 
GISSEL (Bruzual & Chariot 1993) model (1995 version: see 
Chariot, Worthey & Bressan 1996) in which our BCGs form 
in a 1 Gyr burst (with a Salpeter 1955 IMF, over the range 
0.1 < M/Mq < 125) at a redshift Zt = 3, in a Universe with 
Q.M = 1, f^A = 0, and h — 0.5. This can be approximated 
to better than 0.02 mag over the range < z < 1 by a fifth- 
order polynomial in z, given by 



K{z) = -2A76z 



-0.5809z 2 + 14.04z 3 - 



21.10z 4 +9.147z J .(l) 



This is computed by approximating the K band filter by 
a top-hat between 2.0 and 2.45 /»m: the computed K- 
correction changes by ^ 0.02 mag over < z < 1 if in- 
stead we fit the UKIRT K band filter (S.K. Leggett, private 
communication) with a ninth-order polynomial, which accu- 
rately follows variations in its transmission from 1.9 to 2.5 
/im. The model of equation (1) agrees to better than 0.05 
mag with the mean empirical K band K-correction derived 
by Bershady (1995) from the large sample of z ^ 0.3 field 
galaxies constructed by Bershady et al. (1994), which is very 
uniform, having an rms dispersion of ~ 0.05 mag over the 
full range of field galaxy classes. 

To illustrate the insensitivity of K band K-corrections 
to star formation history we show in Fig. j3| the difference 
between the default model of equation (1) and the K(z) 
curves resulting from (a) varying Zf, and (b) varying the 
IMF and form of the burst of star formation. It is clear that 
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Figure 3. The difference between the default K band K- 
correction model of equation (1) and other models, generated 
using Bruzual & Chariot's 1995 GISSEL library (see Chariot 
et al. 1996), all using a Salpeter (1955) IMF and mass range 
0.1 < M/Mq < 125 unless otherwise stated: the oscillating form 
of these AK(z) curves is caused by the use of Chebyshev polyno- 
mials to produce the fit given in equation (1). In (a) we consider 
bursts of star formation of 1 Gyr duration at Z{ = 1 (dotted line) , 
1.25 (dashed line), 1.5 (dot-dashed line), 2 (dot-dot-dot-dashed 
line), 3 (solid lines) and 4 (long-dashed lines): the K(z) curve 
does not change significantly for Zf > 4. In (b) we show bursts at 
Zf = 3 with: (i) 1 Gyr duration (solid line, as in (a)); (ii) expo- 
nential decaying star formation rate with r = 1 Gyr (dotted line) ; 
(iii) instantaneous burst (dashed line); (iv) instantaneous burst 
with mass range 0.1 < M/Mq < 2.5 (long-dashed line); (v) in- 
stantaneous burst with Scalo (1986) IMF (dot-dashed line); and 
(vi) instantaneous burst with Miller & Scalo (1979) IMF (dot- 
dot-dot-dashed line). We assume SIm = lj ^A = and h = 0.5 in 
all cases. The pronounced dips in K(z) at 2=0.7—1.0 for Z{= 1.0— 
1.5 are caused by AGB stars: their appearance is short-lived, and 
pushed to z > 1 if Zf > 1.5. We only have one BCG at z > 0.6, 
so their effect on our results will be negligible. 



the age of the stellar population is the most important factor 
in determining the K-correction, but the variation of the K- 
correction with age only becomes noticable above z ~ 0.4, 
and is pushed to higher redshifts as zt increases: as long as 
the bulk of a galaxy's stars were formed before zt — 1.0, then 
its K-correction is known to 0.05 or better if its redshift is 
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Z < 0.6 (we only have one BCG more distant than that) and 
known to that accuracy out to z ~ 0.8 if Zi > 1.5. Fig. H (b) 
shows that K{z) is very insensitive to the IMF and duration 
of the burst of star formation. 

If we convert the apparent magnitudes in Fig. H to ab- 
solute magnitudes, Mk, using the K-correction model of 
equation (1) then the K-corrected mean absolute magni- 
tude of the full 47 galaxy sample is Mr = —26.40 mag, 
and the rms dispersion about that mean is a = 0.47 mag. 
If we consider the higher-redshift half of the sample, (i.e. 
the 23 galaxies with z > 0.227) then these figures become 
Mk = —26.53 mag and a — 0.30: this reduction in the rms 
dispersion in Mr is significant at the 98.6 per cent level, 
since only 1399 out of 10 5 23-galaxy random subsamples of 
the full BCG sample give a values less than 0.30. Rather 
than assuming that the BCGs were formed at a single red- 
shift, Z{, we could follow Glazebrook et al. (1995) and assume 
that the stellar populations of our BCGs are of a constant 
age, irrespective of the redshift at which they are observed: 
that procedure leads to very similar results as those of the 
K-correction of equation (1) as the age of the BCGs is varied 
over the range ~ 1 — 10 Gyr (Glazebrook et al. 1995). The 
constant zt and constant age assumptions both represent 
crude treatments of the star formation history of our BCGs, 
but the essential point here is that the very weak dependence 
of K(z) to the exact formation redshift or exact age illus- 
trates how insensitive near-infrared K-corrections are to the 
details of the star formation history of galaxies, and, hence, 
show the power of K band observations of early-type galax- 
ies with a narrow intrinsic spread in absolute magnitude to 
act as standard candles. 

It might be objected that the passive evolution mod- 
els discussed in this Section are not appropriate for BCGs, 
since they might have experienced more recent star forma- 
tion, either due to cooling flow activity or merging with other 
cluster galaxies. Allen (1995) studied six cooling flow clus- 
ters, finding large UV/blue continua and strong emission- 
lines, such as Ha and [Oil], which can be interpreted as 
emission from O stars. In order to test for ongoing star for- 
mation we have examined the difference in absolute magni- 
tudes between the EMSS BCGs in our sample which show 
some evidence for star-forming activity in their nucleus and 
the rest. From the EMSS cluster catalogue (Gioia & Lup- 
pino 1994) the BCGs in MS0302.7, MS0419.0, MS0537.1, 
MS0839.8, MS0955.7, MS1004.2, MS1125.3, MS1224.7 and 
MS1455.0 all have [Oil], while MS0015.9 has many galaxies 
with "E+A" type spectra. Despite this, all ten lie very close 
to the mean Hubble relation shown in Fig. ti. A two-sided 
K-S test for BCGs in our sample with and without strong 
emission lines yields a probability of 0.17 that the difference 
in absolute magnitudes between the two populations would 
arise if the two subsamples are drawn from the same parent 
population. This result suggests that the effect of any on- 
going star formation in the centres of BCGs on the K band 
magnitudes integrated over the 50 kpc diameter aperture is 
small. Such an interpretation of the spatial extent of any star 
formation in these systems is consistent with a number of 
studies which have concluded that optical colour gradients, 
indicative of active star formation, are either completely ab- 
sent in BCGs (Andreon et al. 1992) or are confined to their 
inner regions (McNamara & O'Connell 1992). 

Some BCGs are radio galaxies, so another possible con- 
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cern is that some of the K band light from our BCGs comes 
from AGN, rather than from stars (Eales et al. 1997, and 
references therein). Gioia & Luppino (1994) list radio detec- 
tions for only ten of the BCGs in our sample, and only six 
of these are in the high X-ray luminosity clusters in which 
we shall be particularly interested: a K-S test for the ab- 
solute magnitude distributions of our BCGs in high X-ray 
luminosity clusters yields a probability of 0.78 that BCGs 
with and without radio detections are drawn from the same 
parent population. 

To summarise, there is no evidence to suggest we should 
not interpret the K band photometry of our BCGs as arising 
from the passive evolution of a mature stellar population 
formed at z <; 3 similar to that advocated by Ellis et al. 
(1997) and others for cluster ellipticals in general. We shall, 
therefore, adopt equation (1) as our preferred K-correction, 
but will return, in Section 8, to discuss the effect that varying 
this assumed K{z) relation has on the constraints on Om and 
SIa we deduce from the BCG Hubble diagfram. 



5 PROPERTIES OF BCGS AND THEIR HOST 
CLUSTERS 

The scatter in the raw Hubble diagram shown in Fig. g 
should come as no surprise, since it is well known that there 
exist strong correlations between the optical luminosities of 
BCGs and properties both of themselves, and of their host 
clusters, and that it is only after correction for these that 
BCGs are revealed to be good standard candles. In this Sec- 
tion we investigate some of the physical properties of the 
BCGs and their host clusters likely to be correlated with 
BCG luminosity. 

5.1 The BCG structure parameter 

Sandage (1972a,b) found that the luminosities of BCGs cor- 
relate with the richnesses and Bautz- Morgan (Bautz & Mor- 
gan 1970) type of their host clusters, and Hoessel (1980) 
showed how this could be expressed as a correlation between 
L m , the BCG luminosity within a metric aperture of radius 
r m , and a = [dlog(L m )/dlog(r)]| r =r- m , the logarithmic slope 
of the BCG growth curve. Correction for the correlation with 
the structure parameter, a, greatly reduces the dispersion in 
BCG absolute magnitudes, and the L m — a relation forms the 
basis of the recent use of BCGs as streaming velocity probes 
(Lauer & Postman 1994; Postman & Lauer 1995; Hudson & 
Ebeling 1997). 

To measure the structure parameters of our BCGs, we 
computed the growth curve for each BCG from luminosity 
measurements in a set of apertures with radii increasing in 
steps of 1 pixel out to 100 pixels or 30 arcsec. Each curve was 
then modelled using a Hermite polynomial, which was fitted 
over the full range of the growth curve. The value of a was 
then determined for each BCG by evaluating the logarithmic 
derivative of the polynomial fit to its growth curve at an 
angle # m , which corresponds to a radius r m = 25 kpc, and 
is given by 



9m (arcsec) = 0.43 - 



(1 + z) 2 



[1-(1+Z)-V2]' 

for our assumed cosmology, with (Jm = 1, ^a = and h 



(2) 
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Table 1. Seeing disk convolution results: the mean increase in a 
for two values of r m , as a function of redshift. 

r m mean Aa mean Aa mean Aa 

(kpc) 0.2 < z < 0.3 0.3 < z < 0.4 0.4 < z < 0.8 



25 



0.04 
0.02 



0.05 
0.01 



0.10 
0.03 



0.5. There are three principal sources of potential systematic 
error in our a estimates: 

(i) Incorrect background subtraction. In all cases, our 
galaxies are small enough that an adequate portion of un- 
contaminated image was available to estimate the sky level. 
We estimated what effect residual variations in the sky back- 
ground have on the estimated structure parameter by recal- 
culating the sky level using more than one uncontaminated 
part of each image. These tests were carried out on a ran- 
dom subset of 28 BCGs in our sample. The results indicated 
that uncertainties in sky subraction give rise to typical la 
errors of about 0.06 which is ~ 10 per cent on average: this 
dominates the uncertainty in the estimation of the structure 
parameter. 

(ii) The presence of contaminating sources within the 
metric aperture. For 7 sources a star or second galaxy was 
found to lie within the circle corresponding to a projected 
distance of 25 kpc from the centre of the BCG. In these cases 
the contaminating sources were masked out, and the values 
of the pixels in the masked region, (which, in all cases, were 
close to 25 kpc away from the BCG), were replaced with 
the value of the pixel located at the same distance from the 
BCG but diametrically opposite the contaminating source. 
In tests with non-contaminated images, this process intro- 
duced very little additional error. 

(iii) Profile broadening due to seeing. In the inner region 
of the BCG growth curves, the broadening effect of seeing 
causes the a values to be systematically overestimated. Over 
the course of our observations the seeing varied between 0.9 
and 1.4 arcsec, with a median seeing of 1.1 arcsec: a see- 
ing disk of this size is equivalent to a physical size of 5 kpc 
at z — 0.2 and 8 kpc at z = 0.5. To investigate the size 
of this effect we smoothed our images by convolving them 
with a Gaussian filter with a FWHM of 1.1 arcsec and re- 
measured the structure parameters. The mean increase in 
the measured a parameter computed at r m = 8h -1 kpc and 
r m = 25/i -1 kpc in redshift bins is shown in Table 1. It is 
clear from Table 1 that the estimated corrections to the a 
parameter measured in a 50 kpc diameter aperture are small 
compared to other uncertainties, and there is little variation 
with redshift, so it is not necessary to correct our data for 
the effects of seeing. The advantages of using a larger aper- 
ture can also be seen, since the systematic corrections are 
significantly smaller for a 50 kpc aperture compared with 
those for a 16 kpc aperture. The adopted a values for our 
BCGs are tabulated in Table 2. 



5.2 Host cluster X-ray luminosities 

Hudson & Ebeling (1997) have recently argued that the 
residuals in the L m — a relation correlate with the X-ray 
luminosity of the host cluster, so that the L m — a rela- 
tion does not fully remove the environmental dependence 



on BCG luminosity: the potential use of the cluster X-ray 
luminosity and/or temperature in reducing the scatter in the 
BCG Hubble diagram was also discussed by Edge (1991). In 
Table 2 we list the X-ray luminosities of the host clusters 
of our BCG sample. The EMSS X-ray luminosities for the 
45 EMSS clusters are taken from Gioia & Luppino (1994), 
and are quoted for the energy pass band of the Imaging 
Proportional Counter (0.3 — 3.5 keV). We supplement that 
information with the luminosities for two ROSAT clusters 
(R84155 and R843053) which were discovered as part of an 
investigation of large-scale structure in the southern hemi- 
sphere (see Romer et al. 1994). The X-ray luminosities of 
these additional clusters have been corrected for extinction 
by using Hi column densities interpolated from Stark et al. 
(1992), then K-corrected and transformed from the ROSAT 
Position Sensitive Proportional Counter energy pass band 
(0.1 — 2.4 keV) to the EMSS energy range assuming a 6 keV 
thermal bremsstrahlung spectrum: the pass band correction 
isoftheformL x (0.3-3.5keV) = 1.08 x L x (0.1 -2.4 keV). 

5.3 Near neighbour analysis 

To test how well the correlations with a and Lx can ex- 
press the environmental dependences of the BCG absolute 
magnitude, we wanted also to study residual dependences 
on the richness of the host cluster. Since our principal goal 
was to measure BCG magnitudes, we did not image more of 
the host clusters than fell into the field of view of the chip 
while imaging the BCG itself. In the most extreme case, this 
means we can only count near neighbours out to a distance 
of 100 kpc from the BCG: this means that we are looking 
at the central core of the cluster, rather than determining 
the richness of the cluster as a whole, for which one would 
want to count neighbours to some appreciable fraction of an 
Abell (1958) radius. 

We used the STARLINK PISA package (Draper & Eaton 
1996) to detect near neighbours within a projected distance 
of 100 kpc of the BCG centre, by searching for sets of at least 
ten contiguous pixels lying at least 2<r above the sky level 
and used the iraf APPHOT package to estimate their mag- 
nitudes: we removed 'obvious' bright stars and foreground 
galaxies by eye, although few BCGs had one of these falling 
within the counting circle. From the sky noise level in each 
frame we were able to deduce its limiting isophote, which 
we converted into a limiting magnitude. We then computed, 
for each frame, the magnitude difference between the 50 kpc 
BCG magnitude and the limiting magnitude for that frame, 
finding that all frames went at least 4.66 magnitudes deeper 
than the BCG magnitude. Our raw near neighbour count for 
each BCG, N, was then taken to be the number of objects 
lying within a projected distance of 100 kpc of the BCG and 
with a magnitude less than 4.5 magnitudes fainter than the 
50 kpc aperture magnitude of the BCG. 

We then estimated, for each BCG, the number of field 
galaxies expected to fall within the 100 kpc counting circle 
and the 4.5 mag magnitude strip. This was done using a fit to 
a compilation of K band galaxy count data kindly provided 
by J. Gardner (private communication), which comprises 
data from a number of authors: Mobasher et al. (1986), 
Gardner et al. (1993), Glazebrook et al. (1994), Djorgovski 
et al. (1995), McLeod et al. (1995), Metcalfe et al. (1996), 
Huang et al. (1997) and Gardner et al. (1997). These data 
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Table 2. Data for the BCGs in the X-ray cluster sample 
Cluster RA dec z 



Lx 



mBCG 



N 



N 



MS0007.2-3532 


00' l 07 m 14!5 


-35°33'11'.'8 


0.050 


0.517 


11.14 


0.01 


0.49 





-0.50 


MS0015.9+1609 


00 h 15 m 58!3 


+16°09'34'.'0 


0.546 


14.639 


15.58 


0.03 


0.57 


5 


3.15 


MS0301. 7+1516 


03 ,l 01 m 43!3 


+15°15'50?6 


0.083 


0.330 


12.33 


0.08 


0.43 





-1.00 


MS0302.5+1717 


03 ,l 02 m 29!4 


+17°16'47'/6 


0.425 


2.879 


14.65 


0.07 


0.75 


5 


3.89 


MS0302.7+1658 


03' l 02 m 43!2 


+16°58'27'.'0 


0.426 


5.043 


15.01 


0.07 


0.67 


1 


2.51 


MS0419.0-3848 


04 ,l 18 m 59?6 


-38°49'01"2 


0.225 


0.385 


14.67 


0.01 


0.32 





-2.36 


MS0433.9+0957 


04 /l 33 m 58!4 


+09°57'36'.'7 


0.159 


4.335 


13.11 


0.10 


0.67 


1 


3.14 


MS0451. 6-0305 


04 ft 51 m 40 s 5 


-03°05'46'.'0 


0.55 


19.976 


15.74 


0.03 


0.61 


6 


3.94 


MS0537.1-2834 


05 h 37 m 06!8 


-28°34'40 / /6 


0.254 


0.599 


14.97 


0.01 


0.28 


3 


0.43 


MS0821. 5+0337 


08 fe 21 m 33!7 


+03°37'30'.'3 


0.347 


1.328 


14.55 


0.02 


0.76 


1 


2.74 


MS0839.8+2938 


08 ,l 39 m 53!3 


+29°38'16?0 


0.194 


5.348 


13.42 


0.03 


0.65 


7 


6.10 


MS0849.7-0521 


08 h 49 m 4ffl3 


-05°21'36'.'5 


0.192 


1.179 


14.49 


0.02 


0.34 


1 


1.47 


MS0904.5+1651 


09 h 04 m 33?0 


+16°51'15'.'0 


0.073 


0.918 


11.59 


0.03 


0.46 


6 


5.53 


MS0906.5+1110 


09 h 06 m 30!l 


+11°10'39'.'3 


0.180 


5.769 


12.97 


0.03 


0.69 


7 


6.39 


MS0955.7-2635 


09' l 55 m 45?2 


-26°35'56?2 


0.145 


2.039 


14.32 


0.04 


0.13 


6 


2.68 


MS1004.2+1238 


10 h 04 m 12?7 


+12°38'15'.'9 


0.166 


0.925 


13.16 


0.03 


0.81 


1 


3.14 


MS1006.0+1202 


10 h 06 m 07!3 


+12°02'20'.'4 


0.221 


4.819 


13.52 


0.03 


0.74 


9 


8.17 


MS1008.1-1224 


10 ft 08 m 05?4 


-12°25'07?4 


0.301 


4.493 


14.12 


0.04 


0.73 


8 


7.17 


MS1054.4-0321 


10 h 54 m 26?9 


-03°21'25'.'3 


0.823 


9.281 


16.14 


0.03 


0.70 


8 


5.85 


MS1111. 8-3754 


ll h ll m 49!8 


-37°54'55'/8 


0.129 


4.325 


12.16 


0.07 


0.70 


2 


1.61 


MS1125.3+4324 


ll h 25 m 17!9 


+43° 24' 09'.' 5 


0.181 


0.756 


14.94 


0.02 


0.23 


3 


-1.00 


MS1127.7-1418 


ll h 27 m 52!3 


-14°18'19?2 


0.105 


2.786 


12.02 


0.03 


0.90 


6 


5.55 


MS1147.3+1103 


ll h 47 m 18?2 


+11°03'15'.'9 


0.303 


2.304 


14.36 


0.02 


0.61 


3 


1.76 


MS1208.7+3928 


12 ,l 08 m 44!0 


+39°28'19?0 


0.340 


2.030 


14.32 


0.01 


0.55 


3 


1.95 


MS1224.7+2007 


12 h 24 m 42?6 


+20°07'30?0 


0.327 


4.606 


14.61 


0.02 


0.62 


1 


-0.41 


MS1241. 5+1710 


12 h 41 m 31!6 


+17°10'06'.'7 


0.312 


3.411 


14.56 


0.01 


0.44 


5 


3.57 


MS1253.9+0456 


12 h 53 m 54?l 


+04° 56' 25'.' 7 


0.230 


3.143 


13.66 


0.01 


0.56 


4 


3.09 


MS1426.4+0158 


14 h 26 m 26!7 


+01°58'36'.'9 


0.320 


3.707 


14.77 


0.02 


0.58 


1 


2.35 


MS1455.0+2232 


14 h 55 m 00?5 


+22°32'34?7 


0.259 


16.029 


13.91 


0.02 


0.56 


3 


2.01 


MS1512.4+3647 


15 h 12 m 25!9 


+36°47'26?7 


0.372 


4.807 


14.95 


0.01 


0.82 


8 


6.38 


MS1522.0+3003 


15 h 22 m 03?6 


+30° 03' 51? 1 


0.116 


2.347 


12.36 


0.01 


0.53 


7 


6.41 


MS1531. 2+3118 


15 h 31 m 14!l 


+31°18'42'.'2 


0.067 


0.444 


12.07 


0.02 


0.35 


5 


3.95 


MS1532.5+0130 


15 h 32 m 29!4 


+01°30'46'.'3 


0.320 


1.641 


14.80 


0.03 


0.36 


3 


1.31 


MS1546.8+1132 


15 h 46 m 52f0 


+11° 32' 25? 6 


0.226 


2.937 


13.52 


0.03 


0.52 


4 


3.19 


MS1558.5+3321 


15 h 58 m 26?9 


+33°21'40'.'7 


0.088 


1.420 


12.33 


0.02 


0.77 


7 


6.10 


MS1617.1+3237 


16 h 17 m 08?8 


+32°37'52'.'6 


0.274 


1.185 


14.22 


0.01 


0.79 


6 


4.76 


MS1618.9+2552 


16 h 18 m 56?7 


+25°53'22?3 


0.161 


2.241 


13.77 


0.03 


0.94 


7 


5.31 


MS1621. 5+2640 


16 h 21 m 32!2 


+26°41'06?4 


0.426 


4.546 


14.96 


0.03 


0.32 


5 


3.57 


MS2053.7-0449 


20 ft 53 m 44?0 


-04° 49' 24? 7 


0.583 


5.775 


15.96 


0.03 


0.51 


1 


1.70 


MS2124.7-2206 


21 h 24 m 39!4 


-22°07'15?2 


0.113 


1.161 


12.05 


0.03 


0.57 


1 


3.58 


MS2215.7-0404 


22 h 15 m 41!3 


-04°04'25?2 


0.090 


1.196 


12.63 


0.04 


0.60 


1 


2.75 


MS2216.0-0401 


22 h 16 m 04!7 


-04°01'51?9 


0.090 


1.935 


12.23 


0.04 


0.93 


5 


4.23 


MS2255.7+2039 


22 h 55 m 40!7 


+20° 39' 04'.' 2 


0.288 


2.041 


14.37 


0.04 


0.62 


8 


6.67 


MS2301. 3+1506 


23' l 01 m 17!l 


+15°06'49'.'8 


0.247 


3.291 


13.75 


0.04 


0.52 


3 


2.10 


MS2354.4-3502 


23 ,l 54 m 25!9 


-35°02'15?8 


0.046 


0.392 


10.49 


0.01 


0.51 


1 


0.78 


R84155 


01 h 42 m 20!2 


-22°28'39?8 


0.278 


6.97 


13.85 


0.01 


0.55 


5 


4.14 


R843053 


23' l 47 m 47!7 


-24°52'30?9 


0.193 


3.28 


13.73 


0.01 


0.42 





-1.24 



The names, positions, redshifts and X-ray luminosities of the EMSS clusters (with prefix MS) are taken from 
Gioia & Luppino (1994): the RAs and decs, are given in B1950 coordinates and refer to the brightest cluster 
member in the V band; and the X-ray luminosities, Lx, arc in units of 10 44 erg s — 1 . mBCG i s the BCG K 
band magnitude within a 50 kpc diameter aperture, Ak is the K band extinction correction (computed using 
the Afjji maps of Stark et al. 1992, Ry = 3.5 and the mean extinction law of Mathis 1990), a is the structure 
parameter at 50 kpc diameter, while N and iV are, respectively, the raw and corrected BCG near neigh- 
bour numbers. Clusters commonly known by another name are as follows: MS0015. 9+1609 (C10016+16); 
MS0839.8+2938 (Zwicky 1883); MS0904.5+1651 (Abcll 744); MS0906.5+1110 (Abcll 750); MS1006.0+1202 
(Zwicky 2933); MS1127.7-1418 (Abell 1285); MS1253.9+0456 (Zwicky 5587); MS1522. 0+3003 (Abell 2069); 
MS1531.2+3118 (Abcll 2092); MS1558.5+3321 (Abcll 2145); MS1618. 9+2552 (Abcll 2177); MS2255. 7+2039 
(Zwicky 8795); MS2301. 3+1506 (Zwicky 8822); MS2354.4-3502 (Abell 4059); and R84155 (Abell 2938) 
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Figure 4. K band field galaxy count data. The solid line plots the 
fourth-order polynomial fit made to the count data, and used to 
predict the background count for each BCG. The symbols for the 
galaxy counts are as follows: diamonds, Mobasher et al. (1986); 
empty stars, Gardner et al. (1993); filled triangles, Glazebrook et 
al. (1994); crosses, Djorgovski et al. (1995); squares, McLeod et 
al. (1995); circles, Metcalfe et al. (1996); filled stars, Huang et al. 
(1997); and empty triangles, Gardner et al. (1997). 
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Figure 5. The K-z diagram from Fig. fcj showing the location 
of members of quartiles in a: the first quartile is shown by the 
squares, the second by the diamonds, the third by the triangles 
and the fourth by the circles. 



are shown in Fig. W, together with the following fourth order 
polynomial fit to them: 



log 10 «JV(fO» 



oo + a\K + a 2 K + CL3K + a^K , 



(3) 



over the range 11.25 <, K <, 23.5, where (N(K)} is the 
expected galaxy count per mag per degree 2 at magnitude 
K, and the values of the coefficients in eqn (3) are as fol- 
lows: 00 = 26.509, 01 = -8.3590, a 2 = 0.86876, a 3 = 
-3.5824 x 10~ 2 and a 4 = 5.2630 x 10" 4 . This fit has 
X 2 = 128 for 56 degrees of freedom and this is improved 
only marginally by increasing the order of the fit to tenth, 
say. In Table 2 we tabulate, for each BCG, the raw count N 
and the corrected count, N, obtained from it by subtract- 
ing the model background count. In six cases we find that 
the estimated contamination exceeds the raw count N: this 




0.1 



1.0 



Figure 6. The K-z diagram from Fig. ti showing the location 
of members of quartiles in Lx : the first quartile is shown by the 
squares, the second by the diamonds, the third by the triangles 
and the fourth by the circles. 
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Figure 7. The K-z diagram from Fig. fcj showing the location 
of members of quartiles in N: the first quartile is shown by the 
squares, the second by the diamonds, the third by the triangles 
and the fourth by the circles. 

is not unexpected, given the uncertainty in the background 
counts near the limiting magnitudes of the frames (which 
are, typically, 20 Si K ^, 21) and our neglect of clustering in 
the background population, etc. 



6 CORRELATING BCG AND HOST CLUSTER 
PROPERTIES 

The correlations between BCG and host cluster properties 
are graphically illustrated by Figs. H, H andR, which show, on 
the Hubble diagram of Fig. 0, the members of the quartiles 
of, respectively, a, Lx and N. The similarity between these 
three figures shows how all these quantities are broad indi- 
cators of cluster richness, and, in particular, how the fainter 
BCGs at z ^ 0.25 tend to be low a galaxies in poor systems, 
as defined both by the X-ray luminosity, Lx, and by the 
near neighbour number, N. This suggests that there would 
be redundancy in correcting the K — z relation for all three 
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Table 3. Results of Spearman rank correlation analysis. The up- 
per half of the table gives the Spearman rank correlation coef- 
ficient, r^a, while the lower half gives Pab = — log 10 [P(r ab)L 
where P(t"ab) i s the probability that the absolute value of the 
correlation coefficient would be as large as |»"ab| f° r a sample of 
the same size, under the hypothesis that the quantities A and B 
are uncorrelatcd. 





z 


M K 


Lx 


N 


a 


z 


- 


-0.22 


0.61 


0.02 


0.10 


M K 


0.86 


- 


-0.49 


-0.39 


-0.40 


Lx 


5.36 


3.29 


- 


0.34 


0.31 


N 


0.04 


2.15 


1.74 


- 


0.58 


a 


0.28 


2.29 


1.49 


4.73 


- 



quantities, so, in this Section, we perform a non-parametric 
correlation analysis (using the Spearman rank correlation 
coefficient, e.g. Press et al. 1992), to study the relationship 
between BCGs and their host clusters and, hence, to deter- 
mine the best way to correct the BCG Hubble diagram for 
the effects of environment. We consider the five properties 
tabulated in Table 2: (i) redshift, z; (ii) absolute K band 
magnitude, Mk, computed using the K-correction model of 
equation (1); (iii) X-ray luminosity, Lx', (iv) corrected near 
neighbour number, N; and (v) BCG structure parameter, 

Q. 

In Table 3 we present the results of this Spearman rank 
correlation analysis: in the upper half of the table we tab- 
ulate the Spearman rank correlation coefficient, tab, while 
the lower half gives Pab = — log 10 [P(r ab)], where P(r ab) is 
the probability that the absolute value of the correlation co- 
efficient would be as large as |t"ab| for a sample of the same 
size, under the hypothesis that the quantities A and B are 
uncorrelated. The P(rAB) values are computed under the 
assumption that the quantity tabKN — 2)/(l — f"^)] 1 / 2 , 
where N is the number of galaxies in the sample, is dis- 
tributed like Student's t statistic: from Monte Carlo simula- 
tions we judge this approximation to be very good, typically 
estimating Pab to an accuracy of ~0.02. A very similar set 
of results to Table 3 is obtained if we use either the constant 
age K-correction of Glazebrook et al. (1995) or the empirical 
K-correction of Aragon-Salamanca et al. (1993), rather than 
the model fit of equation (1). Note that the use of a rank 
correlation coefficient underestimates the true level of cor- 
relation between Mk and a to the extent that this relation 
is parabolic, as commonly assumed, rather than monotonic 
(see Section 7 below) : none of the other relations tabulated 
in Table 3 are so affected, justifying our use of a rank cor- 
relation analysis to quantify them. 

A number of these correlations are worthy of comment. 
Firstly, the K band absolute magnitudes of our BCGs show 
no evidence of significant evolution with redshift, above the 
passive evolution removed by the K-correction, but we see 
from Table 3 that Mk is strongly correlated with Lx, N 
and a. The strong correlation between z and Lx simply 
reflects the flux-limited nature of EMSS cluster sample, but 
note that the strong correlation between N and a (shown 
in Fig. 0) is not due to near neighbours lying within r m of 
the BCG, since we explicitly excluded such contamination in 
calculating a: instead, it illustrates how well a can express 
the richness of the BCG's local environment, as argued by 
Hoessel (1980). 

These results could, of course, be due to correlations 
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Figure 8. Correlation of N and a. This relation is not simply 
the result of contamination within a radius r m of the BCG, since 
these objects were removed before calculating a. 



with a third quantity: for example, we know that N and a 
are strongly correlated, so the correlation between N and 
Mk, say, could just result from that and a correlation be- 
tween Mk and a. To test that we have computed, for our 
five physical properties, partial Spearman correlation coeffi- 
cients (e.g. Yates et al. 1986) of the form rAB.c, where 



7"AB,C = 



TAB — 7"AC?"BC 



[(!• 



'AC 



)(1 



' BC 



)]l/2 



(4) 



where tab, etc, are the Spearman rank correlation coeffi- 
cients, as before. The significance of the partial rank cor- 
relation coefficients may be calculated using a similar ap- 
proximation as above - that rAB,c[(N — 3)/(l — t"ab,c)] 
is distributed as Student's t statistic - to compute the prob- 
ability, P(rAB,c), of obtaining a partial rank correlation co- 
efficient with absolute value as large as |rAB,c| under the 
null hypothesis that the correlation between A and B re- 
sults solely from correlations between A and C and between 
B and C: Monte Carlo simulations show that this approxi- 
mation works as well for P(r ab,c) as for P(rAB). 

We relegate the full set of partial rank correlation re- 
sults to an Appendix, and note only a few important points 
here. The most striking result is that once the dependence 
on X-ray luminosity has been removed, a becomes indepen- 
dent of z to a high degree: P(r ZQi L x ) = 0.60. We find that 
the N — a relation is not due to separate correlations be- 
tween those two variables and a third: the lowest P(rAB,c) 
for A = N and B = a arises when C — Mk, but, even then, 
there is a probability of less than 4 x 10 -4 that the N — a 
correlation is entirely due to the Mk — N and Mk — a corre- 
lations. This leads to high values of P(r ab,c) when A = N 
and C = a, or vice versa, since these variables are almost 
interchangable. There is some reduction in the significance 
of the Mk — Lx correlation when it is evaluated at constant 
a, but it remains quite strong, consistent with the sugges- 
tion by Hudson & Ebeling (1997) that correction for Lx 
can remove some of the residual environmental dependence 
of BCG magnitudes left after the L m — a correction has 
been made. Finally, it appears that none of the significant 
correlations between pairs of properties arise from separate 
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Figure 9. The L m — a relation for our data set. The filled and 
empty squares denote the BCGs in clusters with X-ray luminosi- 
ties greater than and less than 2.3 X fO 44 erg s — 4 respectively. 
The line shows the best quadratic fit to the data. 



evolutionary trends of the properties with redshift. To sum- 
marise, it is clear that Mk is strongly correlated with both a 
and Lx, and that no further correction for TV need be made 
after the a-correction, due to the strength of the correlation 
between N and a. 



7 THE CORRECTED K-z RELATION 

The relationship between absolute magnitude and a for our 
BCG sample is shown in Fig. H. The relationship diverges 
from a purely linear slope at large values of a and is in 
this respect qualitatively similar to those of Hoessel (1980), 
who used a sample of 108 Abell clusters, and Postman & 
Lauer (1995), who determined a for a sample of 119 BCGs 
at z < 0.05 in the R band. In a similar way to Postman 
& Lauer (1995), we fit the correlation between structure 
parameter and K band luminosity with a parabola, which 
gives a least- squares best fit of 

M K = -23.84- (8.34 ± 1.14)a + (6.14 ± 1.00)a 2 , (5) 

which is also shown in Fig. M The vertical scatter about this 
fit is a = 0.29 mag, and the gradient in the linear part of the 
relation is significantly steeper for our sample than for the 
R band sample of Postman & Lauer (1995). It is clear that 
a parabolic fit is a reasonable description of the data over 
the range 0.1 < a < 1.0, and Fig. M provides the first real 
evidence that the L m — a relation at least flattens at high 
q, rather than continuing linearly with increased scatter, as 
appears the case for the Postman & Lauer (1995) data. 

The filled squares in Fig. M are the BCGs in high-Lx 
clusters. Note that these are all to be found near the peak 
of the parabola, indicating how high-Lx selection picks out 
BCGs with a narrow range of absolute magnitudes, and sug- 
gesting that these magnitudes will be largely unaffected by 
the correcting for the Mk — a relation of equation (5). This 
is graphically demonstrated by Figs. |l0| and y_l|, which show, 
respectively, the relationship between Lx and the raw and 
Q-corrected Mk values. Below host cluster X-ray luminosi- 
ties of Lx — 2.3 x 10 44 erg s _1 , there is a wide dispersion 
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Figure 10. The relationship between X-ray luminosity and raw 
BCG absolute magnitude. The empty squares denote the BCGs 
in clusters with Lx < 2.3 X 10 44 erg s — 1 , while the filled squares 
are the BCGs in clusters above that threshold. Note the much 
wider spread of absolute magnitudes seen for the BCGs in the 
less X-ray-luminous clusters. 
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Figure 11. The relationship between X-ray luminosity and a- 
corrected BCG absolute magnitude. The empty squares denote 
the BCGs in clusters with L x < 2.3 X 10 44 erg s~\ while the filled 
squares are the BCGs in clusters above that threshold. Compar- 
ison of this plot and Fig. id shows how the a correction removes 
most of the scatter in the absolute magnitudes of the low-Lx 
BCGs, while leaving the high-Lx ones virtually unchanged. 



in raw BCG absolute magnitudes, which is largely, but not 
completely, removed by the Mk — a correction, while the 
high-Lx absolute magnitudes are barely changed. In what 
follows we use this threshold of Lx = 2.3 x 10 44 erg s _1 to 
mark the break between a homogeneous population of BCGs 
in high-Lx and a low-Lx population which displays a wider 
scatter in its properties. Intriguingly, the Lx threshold we 
empirically found on the basis of BCG properties is almost 
identical to that Annis (1997) has recently argued marks 
the dividing point in the X-ray properties of EMSS clusters: 
above 2 x 10 44 erg s _1 the X-ray properties are very homo- 
geneous, while below that, Annis detects two populations 
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of clusters, with quite differing properties, one of which he 
argues contains clusters where the X-ray emission is solely 
from the intracluster medium (the only sort of cluster found 
above 2 x 10 44 erg s _1 ) and the other where a significant con- 
tribution comes from a cooling flow. We note that the lack 
of strong correlation between Mk and a for high Lx clusters 
found here, is in contrast to the results found by Hudson & 
Ebeling (1997), who argue that the optical L m — a relation 
actually steepens for high-Lx Abell clusters. However, their 
local Abell sample does not probe luminosities as high as 
the EMSS clusters - only 2 clusters in their sample have 
L x > 2.3 x 10 44 erg s^ 1 in the EMSS passband. 

After correction for the Mk — a relation of equation 
(5) we find that the rms dispersion in absolute magnitudes 
falls from a = 0.47 mag to a — 0.29 magnitudes, while 
the dispersion for the high-Lx subsample (the 26 galaxies 
in clusters with Lx > 2.3 x 10 44 erg s _1 ) about its own 
best-fit Mk — a relation is 0.22 mag, identical to the rms 
dispersion in their raw magnitudes: the best-fit relation for 
the high-Lx sample is M K = -26.18 - (0.87 ± 2.34)a - 
(0.18 ± 1.89)a 2 , indicating the weakness of the dependence 
of Mk on a shown in Fig. U for these BCGs. 

The residuals about the best Mk — a fit for the high-Lx 
sample do not correlate with z, Lx or N, while, for the full 
sample, the only significant correlation is with Lx, similar 
to that reported by Hudson & Ebeling (1997). We also tried 
correcting our raw absolute magnitudes using the same form 
of L m — a — Lx relation recently advocated by Hudson & 
Ebeling (1997) which, translated to the K band, reads 

M K - Co + ciq + c 2 a 2 + c 3 log 10 (L x ) + c 4 Qlog 10 (L x ) 

+ c 5 a 2 log 10 (L x ), (6) 

where Lx is the cluster X-ray luminosity in units of 10 44 
erg s _1 . For the full sample, the best fit values of these 
parameters are Co = —23.26, ci = — 10.66±1.58, C2 = 8.45± 
1.49, c 3 = -3.63±1.13, c 4 = 12.23 ±4.30 and c 5 = - 10.37 ± 
3.97, which means that the terms are in similar proportion 
to those found by Hudson & Ebeling (1997) for the Lauer & 
Postman (1994) optical BCG sample. With this model, the 
rms dispersion of the full sample is decreased slightly to 0.26 
mag, while that for the high-Lx sample actually increases 
slightly to 0.23 mag, as the loss of degrees of freedom arising 
from the use of three more parameters in the fit outweighs 
the marginal decrease in the scatter in absolute magnitudes 
brought about by correction for the residual dependence on 
Lx left after correction for a. 

In principle, galaxy merging in the centre of clusters 
should increase the stellar mass of the BCGs with time. 
Hausman & Ostriker (1978) show that the a parameter is 
a sensitive measure of the merging history of BCGs: there- 
fore, we can determine the effect that the merging process 
has had on our sample by using the a parameter as a tracer 
of the growth of the BCGs through merging. As the Spear- 
man results of Section 6 for the full sample show, once the 
dependences of a on Lx and of Lx on z are removed, a is 
independent of z: there is a 60 per cent probability that the 
correlation between a and z arises from the combination of 
their individual correlations with Lx. This lack of a corre- 
lation between a and z is slightly more strongly confirmed 
for the high-Lx sample of 26 BCGs, as the Spearman test 
for a correlation between a and z for these galaxies gives 
a correlation coefficient with probability of 0.63 of arising 



J- 



from uncorrelated variables. We conclude, therefore, that 
the BCGs in the high-Lx clusters sample form a very homo- 
geneous population, with an intrinsic dispersion in absolute 
magnitude of no more than 0.22 mag, and that there has 
been no significant evolution due to galaxy merging since a 
redshift of z ~ 0.8. It follows that our high-Lx BCG sample 
constitute a set of standard candles ideal for use in classical 
cosmological tests, such as the estimation of Qm and il\. 



8 ESTIMATING ft M AND fi A FROM THE BCG 
HUBBLE DIAGRAM 

The Hubble diagram for a set of standard candles determines 
the relationship between luminosity distance, cZl, and red- 
shift, which may be written (Carroll, Press & Turner 1992; 
Perlmutter et al. 1997) in the general form 

UW\ f*dx [(1 + x) 2 (l + £l M x) - x{2 + x)£l A ]~ 1/2 \ (7) 

where: (i) if Qm + Q,a > 1, then T(x) = sin(a;) and k — 
1 — Am — Qa\ (ii) if Hm + Oa < 1, then T(x) — sinh(a;) and 
K — 1 — Qm — S7a; and (iii) if Qm + &a = 1, then 1(x) — x 
and k — 1. As emphasised by Perlmutter et al. (1997), it is 
important to note that the form of equation (7) means that 
dh depends on J7m and Qa independently, and not solely 
through the deceleration parameter, qo = $7m/2 — Qa - i.e. 
the confidence region in the (Qm, Qa) plane determined from 
the Hubble diagram of a set of standard candles is not par- 
allel to contours of constant qo, except at z <C 1. We prefer, 
therefore, to solve for Q.m and Q.a using our Hubble diagram 
for BCGs in high-Lx clusters, rather than for qo- 

We do that, using the following procedure for each can- 
didate cosmology: 

(i) Select values for Hm and S1a: h = 0.5 is assumed 
throughout. 

(ii) Using equation (7), compute the angle, 6 m , corre- 
sponding to a physical distance r m = 25 kpc at the redshift 
of each of the 47 BCGs: 9 m — t^/cIa, where <La = d L /(l+z) 2 
is the angular diameter distance. 

(iii) Calculate the change in the apparent magnitudes of 
the BCGs, as Attibcg = a(d,A — rfA,ref)/rfA,ref, where the a 
values are those in Table 2, which assume Qm = 1 and 
Qa — 0, and dA,ref is the angular diameter distance in that 
reference model. This first order correction to the apparent 
magnitude should be accurate so long as r m does not vary 
much from the reference aperture with £7m = 1 and Oa = 0, 
and that condition holds true for all physically reasonable 
cosmologies. 

(iv) Compute the ages of the galaxies in the cosmology, 
assuming their formation at a redshift Zi — 3, and compute 
their K-corrections using the Bruzual & Chariot models, 
as described in Section 4. Use the K-corrections for, and 
luminosity distances of, the galaxies in the cosmology to 
derive the BCG absolute magnitudes. 

(v) Perform a quadratic fit to the new Mk — a relation, 
and use it to correct the BCG magnitudes to a fiducial a 
value, which is taken to be the median a. 

(vi) Take the mean of the a-corrected Mk as the estimate 
for Mo, the BCG K band absolute magnitude, and deduce 
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Figure 12. The joint constraints on the value of !!m an d the 
BCG formation redshift Z{ deduced from the Hubble diagram for 
the BCGs in our high-Lx clusters. The dotted line shows the 
locus of the best-fit Cm value for one hundred values of Z{ in the 
range 1 < Z{ < 10, while the solid contours mark the 68, 95 and 99 
per cent confidence contours for Cm and zf. the small amplitude 
random fluctuations in the dotted line are numerical artefacts 
introduced by our method of estimating Cm i while the behaviour 
of all the curves in the range 1.25 < z < 2.25 is caused by the 
transient presence of light from AGB stars in the stellar synthesis 
model used to compute the BCG K-corrections, as discussed in 
Section 4. 



its standard error, from the scatter of the Mk values about 
M . 

(vii) For the redshifts of the 26 BCGs in high-Lx clusters, 
calculate the predicted apparent magnitude, given Mq, and 
the K(z) relation in the cosmology. Using the errors on the 
measured magnitudes of the 26 BCGs in high-Lx clusters, 
and the standard error on Mo, compute \ 2 f° r the compar- 
ison of the 26 BCG apparent magnitudes and the apparent 
magnitudes predicted at their redshifts in the cosmology. 

We report these results for two classes of cosmology, 
with and without a cosmological constant. 



Figure 13. The joint constraints on !!m an d Ca obtainable from 
the Hubble diagram of the 26 BCGs in high-Lx clusters. The 
contours mark the 68, 90 and 95 per cent confidence levels for 
Cm an d Ca jointly! while the the dashed and dotted lines mark 
the models with Ca = and Cm + Ca = 1 respectively. 



narrow intrinsic dispersion in absolute magnitude. Thirdly, 
BCGs are to be preferred over radio galaxies, both because 
the intrinsic dispersion in radio galaxy absolute magnitudes 
appears larger (~ 0.5 mag) and because bias may be in- 
troduced by the significant non-stellar component to the K 
band light (Dunlop & Peacock 1993; Eales et al. 1997). 

If we assume J7a = 0, then we obtain £Im = 0.28 ± 0.24 
if zf = 3. Consistent values are obtained if Zi is varied be- 
tween Zf — 1 (i.e. slightly higher than the redshift of our 
most distant BCG) to Z{ — 10, as shown in Fig. \VA These 
results are consistent with the results of Perlmutter et al. 
(1997), who find Q M = 0.88±g;|g for a tt A = cosmology, 
from the Hubble diagram of the first seven high-z Type la 
supernovae discovered by The Supernova Cosmology Project 
(although our best-fit Qm value is clearly lower than theirs, 
and our confidence interval is much narrower) and with those 
of Carlberg et al. (1996), who find fi M ^ 0.24 ± 0.09 from 
CNOC cluster mass profiles. 



8.1 Friedmann models: £7a = 

Much of the previous work trying to estimate cosmological 
parameters from the Hubble diagrams of distant galaxies 
has assumed J7a = 0, in which case the deceleration pa- 
rameter, go, is given by go — £7m/2. Table 4 summarises 
the go values obtained previously, from a variety of galaxy 
samples, using BCGs and radio galaxies, in both the opti- 
cal and the near-infrared. Several points are worthy of note, 
motivating the use of K band photometry of BCGs in X-ray 
clusters to estimate qo- Firstly, as discussed above, the use 
of K band photometry is to be preferred to the use of optical 
photometry, due to the insensitivity of K band light to star 
formation history of the galaxy. This means that our esti- 
mated qo value is insensitive to uncertainties in the age of the 
stellar population of the BCGs and in the stellar IMF: op- 
tical photometry would produce a greater uncertainty in go 
on both those scores. Secondly, X-ray cluster selection is to 
be preferred over optical cluster selection, since, as we have 
shown, it is very effective at selecting galaxies with a very 



8.2 Friedmann-Lemaitre models: Qa ^ 

Fig. [131 shows the joint constraints on Qm and Qa we ob- 
tain if we allow a non-zero cosmological constant: this is very 
similar to the corresponding plot presented by Perlmutter et 
al. (1997). If we consider only spatially flat models, where 
Om + Oa = 1, we find Clu = 0.55^o'Jg , with a corresponding 
95 per cent upper limit to Qa of 0.73, for zt — 3: as above, 
these figures are very stable, varying marginally over the full 
range z { < 10. Perlmutter et al. (1997) find Q M = 0.94JI° ;^ 
if they assume a spatially flat Universe, which is again con- 
sistent with our results: their results translate into a 95 per 
cent confidence level of Qa < 0.51, slightly tighter than 
ours. From the incidence of gravitational lensing of quasars, 
Kochanek (1996) derived a 95 per cent limit of S1a < 0.66, 
while two recent studies have reported significant detections 
of non-zero Qa values: Henry (1997) finds S7m = 0.55 ±0.17 
from the evolution of the cluster X-ray temperature func- 
tion, and Bender et al. (1997) deduce Qm ~ 0.5, £7a ~ 0.5 
from cluster elliptical Fundamental Plane relations. 
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Table 4. Previous estimates of the deceleration parameter qo from the Hubble diagram of BCGs and radio galaxies. The symbols in the 
waveband column denote whether the measurements were carried out in the optical (O) or near-infrared (1R). Wherever possible the go 
values listed arc those quoted by the authors after correcting for evolution, otherwise the raw values are given. 



Sample 


z range 


Waveband 


Number 


</(> 


Reference 


BCG 


< 0.2 


O 


17 


2.5 ±1.0 


Humason, Mayall & Sandage (1956) 


BCG 


< 0.46 





8 


1.0 ±0.5 


Baum (1957, 1961a,b) 


BCG 


< 0.46 





38 


1.5 ±0.9 


Peach (1970) 


BCG 


0.003 - 0.46 





81 


0.96 ±0.4 


Sandage (1972b) 


BCG 


0.005 - 0.46 





98 


1±1 


Sandage & Hardy (1973) 


BCG 


O.Of - 0.47 





68 


0.33 or -1.27 


Gunn & Okc (1975) 


BCG 


0.04 - 0.46 





67 


1.0 ±0.3 


Sandage, Kristian & Wcstphal (1976) 


BCG 


0.04 - 0.75 





50 


1.6 ±0.4 


Kristian, Sandage & Westphal (1978) 


radio 


0.14-0.95 


1R 


6 


~ 


Lebofsky (1980) 


radio 


0.03 - 1.6 


IR 


83 


0-0.6 


Lilly & Longair (1984) 


radio 


0.016-0.9 


1R 


61 


~0.5 


Lebofsky & Eisenhardt (1986) 


BCG 


0.017-0.147 





116 


-0.55 ±0.45 


Hoessel, Gunn & Thuan (1980) 


BCG 


0.15 - 0.83 





11 


0.3 ±(0.21 -0.39) 


Clowe, Luppino & Gioia (1995) 


X-ray BCG 


0.11 -0.83 


IR 


26 


0.14 ±0.12 


this work 



9 DISCUSSION AND CONCLUSIONS 

This paper has presented the K band Hubble diagram for 
a sample of BCGs in X-ray selected clusters. Our principal 
results are that the brightest X-ray clusters contain BCGs 
with a very narrow spread of absolute magnitudes, and that 
this, coupled with the insensitivity of K band light to the 
details of the star formation processes in, and histories of, 
galaxies, makes samples such as ours ideal for performing 
classical cosmological tests. We find that the BCGs in the 
clusters with high X-ray luminosities (Lx > 2.3 x 10 44 erg 
s _1 in the 0.3 - 3.5 keV passband) have an rms dispersion in 
absolute K band magnitude of no more than 0.22 mag, and 
that leads us to an estimate of Om = 0.28 ± 0.24 if A = 0, 
while, for a spatially flat Universe, we find Qm = 0.55 ±0.15 
and a 95 per cent limit of SIa < 0.73, with only a very weak 
sensitivity to the properties of the stellar population of the 
BCGs over a wide range of redshifts. The stability of these 
results demonstrates that BCGs are excellent standard can- 
dles, competitive with Type la supernovae (which also have 
an intrinsic magnitude dispersion of ~ 0.2 mag: Perlmutter 
ct al. 1997) in the estimation of cosmological parameters Qm 
and Qa'- the further refinement of these two techniques over 
the next few years holds out hope for the determination of 
these fundamental cosmological parameters with high pre- 
cision. 

Finally, the lack of evolution we see in the K band, 
shown by the absence of a correlation between the BCG 
structure parameter, a, and redshift, adds to the growing 
body of evidence that cluster ellipticals are old systems, 
possibly undergoing only quiescent evolution of their stel- 
lar populations (e.g. Dunlop et al. 1996; Bender et al. 1996; 
Ellis et al. 1997). It is intriguing to note that a qualita- 
tively similar lack of evolution is seen in the X-ray luminos- 
ity (Collins et al. 1997, Burke et al. 1997) and temperature 
(Mushotsky & Scharf 1997, Henry 1997) functions. We shall 
return to this question of the relationship between the for- 
mation and evolution of BCGs and their host clusters in a 
forthcoming paper, incorporating recently acquired K band 
imaging data for BCGs in clusters from the deep Southern 
SHARC ROSAT cluster survey (Collins et al. 1997): these 
clusters were selected above a flux limit an order of magni- 
tude fainter than that of the EMSS, enabling us to disentan- 



gle the effects of redshift and X-ray luminosity, in a manner 
impossible with the current data set. 
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APPENDIX A: PARTIAL SPEARMAN TEST 
RESULTS 

In this Appendix we tabulate results from the partial Spear- 
man rank correlation analysis of Section 5 that it would be 
too tedious to include in the main body of the text of this pa- 
per. The partial Spearman rank correlation coefficient rAB,c 
(defined in equation 4 of Section 6) quantifies the strength 
of the correlation between A and B, once the effect of their 
individual correlations with C has been removed. The signif- 
icance of rAB,c may readily be computed. This is the prob- 
ability, Pab.c, that a correlation coefficient with absolute 
value as high as |rAB,c| would arise if the correlation be- 
tween A and B arose solely from their separate correlations 
with C: in common with Yates et al. (1986), we find that 
Monte Carlo simulations indicate that this may be com- 
puted sufficiently accurately using the approximnation that 
tab,c[{N - 3)/(l - r ABjC )] 1/2 is distributed as Student's t 
statistic. 

We tabulate here the rAB,c values for the five quantities 
studied in Section 5: (i) redshift, z; (ii) absolute K band 
magnitude, Mk, computed using the K-correction model of 
equation (1); (iii) X-ray luminosity, Lx; (iv) corrected near 
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Table Al. Results of partial Spearman rank correlation anal- 
ysis for constant rcdshift, i.e. C = z. The upper half of the 
table gives the values of tab C> while the lower half gives 
-Pab = — l°gio[-f( r AB,c)]) where -P(?"ab,c) is the probability that 
the absolute value of the correlation coefficient would be as large 
as |fAB,c| f° r a sample of the same size, under the hypothesis 
that the correlations between A and B arise solely from their 
individual correlations with C 



M K Lx 



N 



Table A4. Results of partial Spearman rank correlation analysis 
for constant corrected near neighbour number, i.e. C = N: format 
as for Tabic Al. 





z 


M K 


Lx 


JV 


a 




- 


-0.23 


0.65 


- 


0.11 


li K 


0.92 


- 


-0.41 


- 


-0.23 


^x 


5.97 


2.33 


- 


- 


0.15 


V 


_ 


_ 


_ 


_ 


_ 



0.31 0.94 



0.49 



Wk 


- 


-0.46 


-0.39 


-0.39 


^x 


- 2.86 


- 


0.42 


0.32 


V 


2.17 


2.45 


- 


0.58 



2.15 1.54 4.67 



Table A2. Results of partial Spearman rank correlation analysis 
for constant absolute magnitude, i.e. C = Mr: format as for Table 

Al. 





z 


M K 


Lx 


JV 


a 


z 


- 


- 


0.60 


-0.76 


0.01 


M K 




- 


- 


- 


- 


Lx 


4.91 


- 


- 


0.19 


0.15 


JV 


0.21 


- 


0.69 


- 


0.50 


a 


0.02 


- 


0.48 


3.45 


- 



Table A5. Results of partial Spearman rank correlation analysis 
for constant structure parameter, i.e. C = a: format as for Table 

Al. 





z 


M K 


Lx 


JV 


z 


- 


-0.20 


0.62 


-0.04 


M K 


0.74 


- 


-0.42 


-0.21 


L x 


5.34 


2.40 


- 


0.21 


N 


0.12 


0.78 


0.79 


- 


a 


- 


- 


- 


- 



neighbour number, N; and (v) BCG structure parameter, 
a. In each table, the upper half lists the partial Spearman 
rank correlation coefficient, rAB.c, while the lower half gives 
Pab.c = — log 10 [P (Vab.c)]- We consider the effect of setting 
C equal to each of these in turn, by which we study the 
significance of correlations between pairs of variables A and 
B after removing the effect of correlations between A and C 
and between B and C. 



Table A3. Results of partial Spearman rank correlation analysis 
for constant X-ray luminosity, i.e. C = Lx'- format as for Table 
Al. 





z 


M K 


^x 


N 


a 




- 


0.12 


- 


-0.26 


-0.13 




0.35 


- 


- 


-0.27 


-0.30 


V 


1.10 


1.15 


_ 


_ 


0.53 


X 


0.40 


1.37 


- 


3.84 


- 
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